tt a tt St Ht tt HH HH HH Ha 
tt a tt a a tt i 


T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 


HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HH HHH 
HHHHHHHHHHHHHHH 
HHHHHHHHHHHHHHH 
HHHHHHHHHHHHHHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 


tt tt tt st 4 Ht Ht Ht Ht Ht Ht HH HH 
AAA a 
| LE) 


LLLLLLLLLLLLLLL 
LLLLLLLLLLLLLLL 
LLLLLLLLLLLLLLL 


Vv ww 
YVwN vVwv 
Ww ww 
ww ww 
Ww ww 
Yow uw 
ww vu 
YVYUVVYVUVUYYVUUWY 
YVYUYVVYUYUUYUUY 
>>>>>>>>>> 
>>> >>>>>>> 
>> 
>> 
>> 
>> 
>> 
>> 
>>>>>>>>>> 
>>>>>>>>>> 
ww 
w et ot te 
=) ne 


—) totes 
ee Ceaioen) 
aooncocaacnnaa 
aoaoncocaaaa 
aa aa 
aa aa 
aa aa 
aa aa 
aa aa 
aa aan 
Aa2QnRbAAAAaaAaaAaaAaaaaa 
aSoaonoanananaaoaaana 
mn Mnnwn nn MMnn 
wn WnMwnn now MMnnw 
wn nw nw nn nw “wn 
mn wn wn Wn ww wn 
mn nn nn nw nw nw 
nn mn wn nw nn “nw 
wn ww nn nw nw wn 
mm nn nw nw nw wn 
YVinnn wn Munn nn 
MMWMnm own Munn wn 
the 
poe 
_— os —e 
ee ee — 
Re ee ee ee ee ag ng ng ag ng Og me Og Og Og 
eee ee te ay ee Od Ey ad On ng ne ne Oe OD Oe 
[od aed ew mom 
me mien) —— 
| ane 
-—e 
8888888888 — 
ay 
a4u 
i) 
a4 
sj4o 
au 
ju 
ssssssssss ee 
ad ad od hd I od od I I I I I 


YL 
> 
aot 
a 
” 
= 
oO 
nal 
& 
f=) 
ao 
- 
aw 
=) 
om 
a 
o 
oe 


FS 
- COMPLEX / COMPLEX DIVISION ROUTINE 


; Detailed Current Edit History 
PLEX giving COMPLEX result 
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Vitte grssp ive - COMPLEX / COMPLEX DIVISION ROUTINE 
IDENT /1-005/ ; File: OTSDIVC.MAR 


SRAAAACHSKAAAAAHTAAAKAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAARAEAAEAAAAEAEAAAAAAETEAAEEAEe 


COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


® 
te 
® 
® 
e 
® 
S SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED * 
Y IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE * 
LUSION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER * 
IES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY * 
ER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY * 
* 
® 
. J 
® 
® 
® 
® 
® 
® 
® 
® 


® 
* 
® 
® 
ex 
oa 
cd 
& 
w 
* TRANSFERRED. 
® 
® 
® 
ee 
te 
& 
® 
® 
® 
® 


Ssz2= 


THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
orppokat ite NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 


LARA ARARAAARASALARAAAALELALEE SALA RSE LESS R SEAS RE EARS ERAS ARSE RSE SASS RRR SAAS SO 


FACILITY: MATH LIBRARY 
; + 
3; ABSTRACT: 
Perform COMPLEX division 
VERSION: 0 
HISTORY: 


AUTHOR: 
Jonathan Taylor, 28-Jun-77: Version 0 


MODIFIED BY: 


5 
gTsspive = COMP LEX / COMPLEX DIVISION ROUTINE 16-SEP=1 38+ M1: 3: 4 ee Macro V04-00 
HISTORY ; Detailed Current Edit History 6-SEP-19 MTHRTL.SRCJOTSDIVC.MAR; 1 


09 4 -SBTTL HISTORY : Detailed Current Edit History 
1 

; Edit History for Version 0 of OTSDIVC 

4 oe 

5 ; 0-04 = Change FORSFLAG_JACKET to MTHSFLAG_JACKET. TNH Bs 16-NoU=78 
000 § ; 1-001 - Update version number and copyright novice. JBS_16-NOV-78 
000 3 e's - lee Late MTHJACKET at assembly time. “9 ¢=7 
0000 8 ; 1-003 = Add “’."’ to the PSECT directive. JBS IS ece 7 
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5 
COMPLEX DIVISION routine 


-SBTTL DECLARATIONS 


: INCLUDE FILES: 
: MTHJACKET.MAR 


: EXTERNAL SYMBOLS: 


; none 
; MACROS: 
: none 


t PSECT DECLARATIONS: 
7 -PSECT _OTSSCODE 


: EQUATED SYMBOLS: 


acce 
or 
or 


OWN STORAGE: 
none 


1 mii . 91:36:36 (tg Macro V04-00 


1 


MTHRTL. SRCJOTSDIVC.MAR; 1 


; Math jacketing macro 


PIC, SHR, 


LONG, EXE, NOWRT 
Program section for OTS$ code 


real part of dividend 
— part of dividend 
real part of divisor 
imag part of divisor 


OTS 
Tat 


5 
LEX DIVISION ROUTINE” 16-SEP-1964 01:52:56 YAK/VIS Macro VO4-00 


-SBTTL COMPLEX / COMPLEX giving COMPLEX result 


a COMPLEX / ¢ 


Etat 
S 2&3 
wm 


ve 
FUNCTIONAL DESCRIPTION: 
OTSSDIVC = COMPLEX / COMPLEX giving COMPLEX result 


The COMPLEX result is computed as follows: 


1) Let (a, b) represent the COMPLEX dividend. 
¢ Let (c, d) represent the COMPLEX divisor. 
>) Let (r, 1) represent the COMPLEX quotient. 


SSSARITH {fs Hb err by hardware on floating zero divide 


SSSARITH GNALed by hardware on floating overflow 


185 § 
i 
1 7 ; 
106 ; 
: 5 3 
§ ; 
i 
0 Hf : Then: 
6 118 : r= (ac + bd) / (cc + dd) 
iW $ i = (be = ad) / (cc + dd) 
000 116 : CALLING SEQUENCE: 
99 118 3 Complex_quotient.wfc.w = OTS$DIVC(dividend.rfc.v, divisor.rfc.v) 
208 0 : INPUT PARAMETERS: 
00 1 § : Dividend and divisor parameters are represented as standard 
O60 : ? : FORTRAN COMPLEX numbers and are CALL BY VALUE. 
9000 125 : IMPLICIT INPUTS: 
800 «159; — 
000 1 8 t OUTPUT PARAMETERS: 
900 9: NONE 
000 131 : IMPLICIT OUTPUTS: 
SES len 
99 4 : FUNCTIONAL VALUE: 
09 6 : The COMPLEX value returned is (a, b) / (c, d) 
000 1 2 + SIDE EFFECTS: 
0 1 ; 
° 140 ; 
141; 
00 142; 
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K 
COMPLEX DIVISION ROUTINE at Tae 9 91:36:36 ee Macro Vv04-00 
OMPLEX giving COMPLEX result -SEP-1984 11:27: MTHRTL.SRCJOTSDIVC.MAR; 1 
144 -ENTRY OTSSDIVC, “M<R2,R3,R4,R5> ; disable integer overflow 
145 MTHSFLAG_JACKET ; establish math error handler 


MOVAB G*MTHSSJACKET_HND, (FP) 
; set handler address to jacket 


3; handier 
148 MOVQ CCAP), RS 3; R4/RS = divisor 
108 3; Perform scaling of all operands before division 
1 : : EXTZ2V #7, #8, R4, Re : R? = R4<exp> 0,1,...377 
1 EXTZV #7, #8, RS, R 3; R1 = RS<exp> 0,1,...377 
1 § CMPW R1, R ; R1 = MAX (c<exp>, d<exp>) 
1 BGTR 2s 
154 MOV R2, R1 
155 2$: MNEGB R1, R1 ; R1 = scaling exponent 0,377,376,... 
138 ROTL #7, RI, RI ; build a floating scale factor 
1 ; scale all operands 
158 MULF R1, R4 s¢ 
159 MULF R1, R5 : d 
160 MULF3 R11, aC(AP), R2 ia 
19] MULF b(AP), R1 3 6 
188 MULF R2,R4,R0 ; RO = ac 
164 MULF3 R1, R5, RB : R3 = bd 
165 ADDF Rie RO 3; RO = ac+bd 
166 MULF R5, R2 3; R2 = ad 
16 MULF R4, R1 : R1 = be 
168 SUBF R2, R1 ; R1 = be - ad 
169 MULF R4, RS ; RO = cc 
170 MULF R5, R5 ; RS = dd 
171 ADDF R4, RS 3; R5 = cc + dd 
v6 DIVF RS, RO 3; RO = (ac+bd) / (cc+dd) 
4, DIVF R5, R1 3; R1 = (be-ad) / (cc+dd) 
175 RET s (RO, R1) = (r, 1) 
178 
17 END 


5 
OTSSOIVC = COMPLEX / COMPLEX DIVISION ROUTINE 
Symbol table 


A z 
C = i 
MTIHSSJACKET_HND eteeeeee x 4 
OTSS$DIVC 00000000 RG 


bower oa wen aos ene sen h 


: Psect synopsis ! 


Tee em ee mmo e eae $+ 


PSECT neme Allocation PSECT No. Attributes 


- ABS, 00000000 ( 0.) 00 ¢ QO.) NOPIC USR CON ABS LC 
-OTSSCODE 00000059 ¢ 89.) O01 ¢ 1.) PIC USR CON REL LC 


teow eee een eeeee Sewer esas eee t 


! Performance indicators H 


gene nena r seco ee ene een ene = 


L NOSHR NOEXE NORD NOWRT NOVEC BYTE 
L SHR EXE RD NOWRT NOVEC LONG 


Phase Page faults CPU Time Elapsed Time 


Initialization 38 00:00:00.09 0:00:01.09 
Command processing 123 0:00:00.63 0:00:04.95 
Pass 82 0:00:00.65 Bp 88! Ba “Ba 
Symbol table sort 0 B90: 80 +80 00:00:00.00 
Pass 2 46 bs Sa 8 5 00:00:01.97 
Symbol table output § 00:00:00.02 00:00:00.02 
Psect synopsis output 00:00:00.01 00:00:00.01 
Cross-reference output 0 ab Bes a 00:00:00.00 
Assembler run totals 286 00:00:01.97 00:00:11.69 


The working set Limit was 750 pages. 

3048 bytes (6 pages) of virtual memory were used to buffer the intermediate code. 

There were 10 pages of symbol table space allocated to hold 6 non-local and 1 local symbols. 
237 source Lines were read in Pass 1, producing 11 object records in Pass 2. 

1 page of virtual memory was used to define 1 macro. 


ese eter ese me sen ewe ecm een ae } 


Macro Library name Macros defined 


0 GETS were required to define 0 macros. 
There were no errors, warnings or information messages. 


"$$EP=1986 91:25:39 EAthntc. seedorspivesmar;1 "#8° 8) 


MACRO/ENABLE =SUPPRESSION/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:OTSDIVC/OBJ=OBJ$:O0TSDIVC MSRCS$:MTHJACKET/UPDATE=(ENHS:MTHJACKET) MSRCS: | 
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